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At the base of all electron-atom scattering and ionization problems is the 
fundamental, yet unsolved, three-body problem of an electron interacting with atomic 
hydrogen. This problem occupies a special place in the set of unsolved problems of 
interest to physicists due to its fundamental nature in the realm of atomic physics. It 
represents a class of Coulomb three-body problems which includes electron interaction 
with the single positive ion of helium, and hence the problem of helium double 
photoionization. 

For heavier atoms the complexity of the Coulomb three-body problem may be 
masked by the collective behaviour of the many target electrons. Similarly, for high 
incident electron energies the complicated role played by the long-ranged Coulomb 
interaction is also somewhat hidden. The problem exhibits all of its complexities at 
energies a little above the ionization threshold for the simplest atomic target, namely 
hydrogen. Here we have the possibility of exciting a countably infinite number of the 
hydrogen discrete states as well as the three-body continuum of two very slow strongly 
interacting electrons. In this Letter we consider the e-H problem at the incident electron 
energy of 15.6 eV, i.e. only 2 eV above the ionization threshold. 

To solve the e-H problem at a total energy E (presently 2 eV) and spin S — 0, 1 
means to correctly predict all of the possible scattering amplitudes /^(k) for discrete 
excitation of target eigenstates with energy e n i < with e n \ + k 2 /2 = E, and ionization 
amplitudes / 5 (kA,k B ) with k\/2 + k 2 B /2 = E. For the discrete transitions the close- 
coupling methods have proved to be the most successful, particularly at low energies. 
These rely on expanding the total wave function in a set of orthonormal states. From the 
landmark work of Yamani and Reinhardt (1975), followed by Broad (1978), Stelbovics 
(1989) and others, it became clear that the set of orthonormal states obtained by 
diagonalising the target Hamiltonian in a Laguerre basis formed an unusual equivalent- 
quadrature rule. Thus-obtained states provide a quadrature rule that incorporates both 
the infinite set of true target discrete states and the true target continuum. This is 
an immensely powerful result and forms the basis of the convergent close-coupling 
(CCC) method for the calculation of electron-atom scattering (Bray and Stelbovics 
1992, Bray 1994, Fursa and Bray 1995). The idea relies on simply increasing the number 
of expansion states N = Z); max Ni until convergence in the amplitude of interest f^ t N (k) is 
obtained to an acceptable accuracy, just like with standard numerical quadrature. This 
approach has proved very successful for the discrete transitions at all energies. In the rare 
case of substantial discrepancy with experiment (Bray and Stelbovics 1992) subsequent 
new measurements were found to be in agreement with the CCC theory (Yalim et al 
1997). 

Obtaining reliable scattering amplitudes for the discrete transitions is a good start, 
but what about ionization? The square-integrable expansion-states 0^ (/ < l max ,n = 
1, . . . , Ni), obtained by diagonalising the target Hamiltonian in a Laguerre basis of size 
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N h have both negative and positive energies e^. With increasing N the negative-energy 
states converge to the true eigenstates (e^ — > e n i, 0^ — > (p n i), and the positive-energy 
states yield an increasingly dense discretization of the continuum. By summing the 
integrated cross sections, obtained upon solution of the close-coupling equations, for 
just the positive-energy states yields excellent agreement with the measurements of the 
e-H total ionization cross section (TICS) (Bray and Stelbovics 1993, Kato and Watanabe 
1995, Scott et al 1997). Though this is the least informative ionization process it is an 
encouraging first step. The question is: do the scattering amplitudes for the excitation 
of the positive-energy 0^ contain all of the detailed ionization information? 

Before proceeding further let us define some convenient notation. Suppose we are 
interested in describing an experiment where the two outgoing electrons have momenta 
\la and k# with k B < k A - When performing the diagonalizations we ensure, by varying 
the exponential fall-off parameter A; (Bray and Stelbovics 1992), that for each / we have 
a state <pn B i whose energy is e^ B i = k B /2. We will refer to these states collectively as 
4> B with energy e B = k 2 B /2. Though it is rarely practical, let us further suppose that 
the same diagonalizations have resulted in states 0^ whose energies e^ A i = k\/2. We 
shall collectively refer to these states as 0^ with energy = k\/2. Similarly, the 
scattering amplitudes f^ N (k) for the excitation of the states <fi B and 0^, arising upon 
solution of the TV-state close-coupling equations, we write as f B N (k A ) and f A N (k B ), 
respectively. Note that the close-coupling formalism ensures that the total wave function 
is expanded explicitly antisymmetrically using the states 0^, with the arising equations 
solved separately for each total spin S. Thus, each /^(k) may always be thought of as 
a combination of the direct F and exchange G amplitudes, eg. /^j (k) = F + (— 1) S G 
for hydrogen. The close-coupling boundary conditions assume that only one electron is 
ever allowed to escape to true infinity, asymptotically as a plane wave. It is helpful to 
keep in mind that the energies e B < e A are symmetrically on either side of E/2, and 
that the summation of the integrated cross sections to obtain TICS includes both sets of 
amplitudes f B N (k.A) and f^ N (k B ) combined as cross sections. For equal-energy sharing 
e B = e A = E/2, which we consider as the limit e B — > 

The work of Bray and Fursa (1996 a) attempted to provide a correct interpretation of 
the already calculated positive-energy-state scattering amplitudes, with some surprising 
and controversial results. It was shown that the (e,2e) ionization amplitudes may be 
defined from the f§ N (k.A) by 



where k^ is a Coulomb wave (in the case of H target) of energy k B /2 = e B . This 
definition is in fact a simplification of the pioneering work of Curran and Walters 
(1987). The overlap has the effect of changing the unity normalization of <p B to that 
of the true continuum, as well as introducing a one-electron Coulomb phase. The 
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controversy (Bencze and Chandler 1999) arises not from the above definition, but from 



The second term above looks like an exchange term, but it is not. The amplitudes 
f§ N (k.A), and hence f SN {kA, kg) are already a coherent combination of their own direct 
and exchange amplitudes as determined by S. The two terms have very different origin. 
The amplitudes /^(k^) arise from the excitation of the states with the boundary 
condition that the "k^" electron exits as a plane wave totally shielded from the ion by the 
bound 0^ electron. For < k A /2 this is the physically sound shielding approximation, 
as used in the Born approximation where the slow electron is modeled by a Coulomb 
wave and the fast one by the plane wave. However, the boundary conditions for the 
amplitude /^(kg) are unphysical (low-energy outgoing plane wave shielded by a higher 
energy bound state). Yet, these two theoretically distinguishable amplitudes correspond 
to the same ionization process since E = e A + e^. 

From (0) we see that close-coupling yields twice as many amplitudes as we may 
expect from formal ionization theory. In the often used language of direct and exchange 
amplitudes we have two such pairs f SN (\s.A, kg) = Fi + (— l) s Gi and / 5Ar (kg, k^) = 
F 2 + (— 1) S G2, which are very different for ^ Note, there is no symmetrization 
relation between the close-coupling theory calculated /^(k^, kg) and / SiV (kg, k^) as 
claimed by Bencze and Chandler (1999). In forming the TDCS we have FiFi, GiGi 
and cross terms G^F^ generally very different for each i = 1,2. A careful numerical 
study of the problem led to the suggestion that with increasing N the second term 
in (0) and hence both F 2 and G 2 converge to zero (Bray 1997). This allows for 
consistency with formal ionization theory except that the f SN (k.A, kg) are obtained 
only for e# < e^. However, for finite N a consistent interpretation (compatible with 
the definition of TICS) of the close-coupling approach to ionization requires the use 
of both terms. A further consequence of the numerical study (Bray 1997) is that the 
close-coupling method is unable to obtain convergence to a satisfactory accuracy in 
the singly differential cross section (SDCS) whenever the true SDCS at equal energy- 
sharing is substantial. Nevertheless, it was argued, that if the true SDCS was known 
then accurate angle-differential ionization cross sections could still be predicted. Here 
we test this claim at just 2 eV above threshold, where the SDCS may be reasonably 
assumed to be approximately flat (Roder et al 1997a). 

The concept of convergence with increasing iV = Y^i<i max Ni involves both the 
increase of Z max and iVj. We performed a series of calculations for various N. The 
ones presented may be conveniently denoted by CCC(iVo, / max ) with Ni = N — I. To 
examine the rate of convergence we present two vastly different calculations CCC(20,5) 
and CCC(13,4), which require approximately 2Gb and 500Mb of computer RAM, 



the subsequent use of (|l|) to define the triply differential cross section (TDCS) by 




^ = | / ^(k A ,k B )| 2 + | / ^(k B ,k^)| 2 . 
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respectively. In both cases the Laguerre exponential fall-off parameter was A; ~ 0.6 
with the variation performed to ensure that for each / there was a state 4>n i with 
energy e„ t = e„ t =l eV. In the present equal energy-sharing case the two terms in (fj) 
are evaluated using the same set of amplitudes, assuming a continuous limit of e B — > e A . 

The first test of the calculations is the comparison of the total ionization cross 
sections (TICS) and its spin asymmetry Aj with the highly accurate measurement (Shah 
et al 1987) of TICS 1.08 (10" 17 cm 2 ) and the A 7 « 0.5 measurements (Fletcher et al 
1985, Crowe et al 1990). The CCC(20,5) and CCC(13,4) results for the TICS, Aj are 
1.18, 0.50 and 0.91, 0.51, respectively. Thus, we see that both calculations attribute 
approximately the correct amount of electron flux to the two spin ionization channels. 
The TICS results from other calculations typically varied around the experimental value. 
The reliability of various close-coupling based theories for the calculation of the TICS 
at low energies has been discussed in detail by Scott et al (1997). The difficulty of the 
problem of obtaining accurate ionization amplitudes at this energy is indicated by the 
fact that the total cross section is more than forty times bigger, 4.7, 4.8 and 4.6±0.1 
(10- 16 cm 2 ) respectively for the CCC(20,5), CCC(13,4) and experiment of Zhou et al 
(1997). 

Next we consider the energy distribution within the ionization channels, i.e. the 
SDCS, defined by 

da SN r 

— (e) = y dn A dn B \f N (k A ,k B )\ 2 . (3) 

The TICS erf is obtained by performing the integration 

^=j[*|(«) (4) 

The integral in (4) is equivalent to the sum of the integrated cross sections for the 
excitation of the positive-energy states. The step function hypothesis (Bray 1997) says 
that the second term in (5) converges to zero with increasing N . The origin of the 
two terms in (^) are the two terms in (5). We think of the second term as numerical 
"left-overs" from an incomplete convergence with N, due to its minor contribution (past 
1 eV) to the TICS. Integration of (|2|), for a given secondary energy e, over the angular 
variables yields the integrand of (5). 

In figure [I] the spin-averaged SDCS are presented. We see that there is no 
convergence in the CCC(20,5) and CCC(13,4) results, though the integral of both is 
much the same. The step function CCC(oo, 5) is an estimate of what the CCC-calculated 
SDCS would converge to for Ni — > oo (there are no problems in obtaining convergence 
with increasing / max ). In other words, we assume that at this low energy the true SDCS 
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Figure 1. The singly differential cross sections arising in the CCC(iVo, i max ) (see 
text) calculations. The step function labeled by CCC(oo,5) is an integral preserving 
estimate. 

is approximately flat. Since the close-coupling theory is unitary we cannot have double 
counting of the TICS, and hence suppose that with increasing iV the SDCS defined in (|3|) 
becomes non-zero only for < e < E/2. In experiment the observed SDCS is symmetric 
about E/2 with the TICS being obtained upon integration to E/2. Comparison with the 
experimental SDCS requires both terms of (5). For the substantially asymmetric energy- 
sharing kinematics only the first term contributes significantly, but both are necessary at 
equal energy-sharing. From figure [I] it is clear that the angular distributions determined 
by (0) will be much too small in magnitude. In order that the integration of (0) over 
the angular variables, the endpoint of the integrand in (5), yielded the estimated SDCS 
of 1.08 (10~ 17 cm 2 /eV) we will multiply the equal energy-sharing CCC-calculated TDCS 
by 1.08/(0.2 x 2) = 2.7. 

In figure ^ we present the TDCS calculated by the two CCC models and compare 
these with experiment and the previously overall best agreement-yielding theory, the 
distorted- wave Born approximation (DWBA) of Jones et al (1992). The relative 
measurements were initially presented by Brauner et al (1991), but were remeasured and 
put on the absolute scale, with an estimated 35% uncertainty, by Roder et al (19976). 
The DWBA calculations (Jones et al 1992) work relatively well at this low energy since 
they utilize the effective charge formalism of Rudge (1968) in the distorting potentials. 
For an example of a more common DWBA approach and the 3C theory see Rouet et al 
(1996) and Brauner et al (1991), respectively. 

In the TDCS figure we use the convenient, for the coplanar geometry, convention 
that the negative scattering angles are on the opposite side of the incident beam [z- 
axis). For best visual comparison with the rescaled CCC calculations we have multiplied 
all of the experimental values by the single constant of 0.45. Having done so, we 
see excellent agreement between the two CCC calculations and experiment for all 
geometries, which is of considerable improvement on the comparison with the DWBA 
calculation. The quality of the agreement gives us confidence that the rescaling of the 
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Figure 2. The coplanar triply differential cross sections, in the indicated geometries, 
for electron-impact ionization of atomic hydrogen with 1 eV outgoing electrons. 
Absolute experiment of Rodcr et al (19976) has been scaled by a factor of 0.45 for 
best visual fit to the rescaled CCC data, see text. The DWBA calculations are due to 
Jones et al (1992). 
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experiment has brought it into consistency with the estimated SDCS value at 1 eV of 
1.08 (lCr 17 cm 2 /eV). Should the true SDCS prove to be a little convex(concave) then 
the experimental rescaling should be done by a factor a little greater(smaller) than 0.45. 
Perhaps the experimentally determined normalization is an indication that the SDCS is 
more convex than concave. As a consequence, we do not believe that the theory of Pan 
and Starace as presented by Roder et al (19976) is a factor of two too low at 15.6 eV, 
and may indeed be accurate at all energies. Though not presented they are almost 
indistinguishable from the 6 AB = 180° rescaled CCC(20,5) TDCS. 

Let us turn specifically to the case where the two detectors are kept 9ab = 80° 
apart. Though no experiment is available for this case we present it because it shows 
the greatest difference between the two CCC calculations, but is still experimentally 
measurable. In fact, smaller Qab geometries yield even greater differences. Such 
geometries, first suggested by Whelan et al (1993), are an excellent test of the CCC 
theory because the cross sections fall rapidly with decreasing 6ab- We see that the 
bigger calculation yields the smaller cross section for Bab — 80°. This is an important 
indication of how well the CCC theory is working. For the other presented cases the fact 
that the shapes of the two calculations are much the same, even though one requires four 
times as much computational resources as the other, suggests rapid shape convergence 
for the largest cross sections. On the other hand, almost identical overall magnitude 
suggests that convergence to the true correct SDCS is extremely slow. 

So how is it that the CCC theory yields such good TDCS angular distributions? To 
help answer this question let us have a look in more detail at the symmetric geometry. 
Given the good agreement between CCC(13,4) and CCC(20,5) TDCS one would imagine 
that one may readily interchange the partial wave amplitudes of ([I]) (kl\<f>^) f^ N (k) in 
the two calculations. The curve labeled by CCC (mix) was generated by taking the 1 eV 
I = 1 partial wave amplitude of the CCC(20,5) calculation and using it with the other I 
CCC(13,4) amplitudes. Whereas one may reasonably expect the CCC(mix) calculated 
TDCS to be between the other two, it differs substantially when the two electrons emerge 
close together. This is an indication of the importance of treating all partial waves in 
a consistent manner. The Laguerre basis choice Ni = No — I with similar A/ results in 
much the same integration rule over the true continuum for each I. In other words, the 
number of positive energy states and their separation is similar for each I. We also use 
the same set of states for each partial wave of total orbital angular momentum J. Thus, 
for each J, the error in the energy distribution is also very similar for each I, and this 
is why the CCC(iVo, / ma x) calculations yield good TDCS angular distributions whose 
magnitude is in error by a single constant. 

What have we learned from this and preceding studies? The CCC approach to e-H 
scattering has not fully solved this Coulomb three-body problem. Given the complexity 
of the problem it is not surprising that the close-coupling approach should run into an 
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intractable problem. Whereas we are confident of obtaining accurate discrete scattering 
amplitudes ab initio at all energies, not so for the ionization amplitudes. Accurate 
ionization amplitudes may require too many states, depending on the incident energy, 
for practical implementation of the CCC theory However, we have suggested two 
empirical prescriptions that still allow for the CCC-calculated ionization amplitudes to 
be useful and predictive, though with some uncertainty. The first, demonstrated here, 
ensures rapid convergence in the angular distributions. It relies on taking a similar 
quadrature rule in the continuum for all target-space /, total orbital angular momentum 
J and total spin S. Defining Ni — N — I with A; « A for each J and S achieves 
this. There may be other more efficient approaches that use a different basis as in say 
the intermediate-energy R-matrix method (Burke et al 1987). A sensible choice for A 
is also important. The second prescription, necessary at low energies when the true 
SDCS is substantially large at E/2, is that of rescaling the cross sections according 
to the ratio of the estimated true SDCS and the close-coupling-calculated SDCS. At 
high-enough energies no such rescaling is necessary and most aspects of the problem 
may be obtained accurately fully ab initio (Bray and Fursa 1996 a, Bray and Fursa 
19966). Given the general structure of the true SDCS at low energies, and that the close- 
coupling based theories obtain the correct TICS, estimating the true SDCS is likely to 
yield only a minor error. Furthermore, in some cases accurate SDCS are available from 
experiment and other theories. With these two empirical prescriptions the close-coupling 
approach to ionization has practical application at all incident energies, energy-sharing, 
and geometries of the two detectors. 

There are a number of opinions relating to improving the close-coupling-calculated 
ionization amplitudes in an ab initio manner such as by matching the calculated total 
wave function to the correct asymptotic three-body boundary conditions. However, we 
note that once the close-coupling equations have been solved the electron flux has been 
incorrectly distributed within the ionization channels. This information will be hidden 
in the total wave function and we suspect will require empirical correction of some kind 
prior to matching. 
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